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ABSTRACT 

Microlensing is increasingly gaining recognition as a powerful method for 
the detection and characterization of extra-solar planetary systems. Naively, 
one might expect that the probability of detecting the influence of more than 
one planet on any single microlensing light curve would be small. Recently, 
however, Griest & Safizadeh (1998) have shown that, for a subset of events, 
those with minimum impact parameter u m j n ~0.1 (high magnification events), 
the detection probability is nearly 100% for Jovian mass planets with projected 
separations in the range 0.6-1.6 of the primary Einstein ring radius Re, 
and remains substantial outside this zone. In this Letter, we point out that 
this result implies that, regardless of orientation, all Jovian mass planets 
with separations near 0.6-1.6-Re dramatically affect the central region of the 
magnification pattern, and thus have a significant probability of being detected 
(or ruled out) in high magnification events. The probability, averaged over 
all orbital phases and inclination angles, of two planets having projected 
separations within 0.6-1.6Re is substantial: 1-15% for two planets with the 
intrinsic orbital separations of Jupiter and Saturn orbiting around O.3-1.OM0 
parent stars. We illustrate by example the complicated magnification patterns 
and light curves that can result when two planets are present, and discuss 
possible implications of our result on detection efficiencies and the ability to 
discriminate between multiple and single planets in high magnification events. 



Subject headings: gravitational lensing, planetary systems 



submitted to the Astrophysical Journal Letters: March 24, 1998 



- 2 - 



1. Introduction 

A planetary microlensing event occurs whenever the presence of a planet creates 
a perturbation to the standard microlensing event light curve. These perturbations 
typically have magnitudes of ~ 20% and durations of a few days or less. First suggested 
by Mao & Paczyhski (1991) as a method to detect extra-solar planetary systems, the 
possibility was explored further by Gould & Loeb (1992), who found that roughly 15% 
of microlensing light curves should show evidence of planetary deviations if all primary 
lenses have Jupiter-mass planets with orbital separations comparable to that of Jupiter. 
Although these probabilities are relatively high, the use of microlensing to discover planets 
was largely ignored since in order to detect the primary events the microlensing survey 
teams must monitor millions of stars in very crowded fields, resulting in temporal sampling 
that is too low (~ lday) and photometric errors that are too high (~ 5%) to detect most 
secondary planetary deviations (Alcock et al. 1997a). 

Recently, the situation has changed dramatically as the real-time reduction of the 
survey teams has enabled them to issue electronic "alerts," notification of on-going events 
detected before the peak magnification (Udalski et al. 1994, Pratt et al. 1996), allowing 
other collaborations to perform special purpose observations of the alerted events. These 
additional observations include denser photometric sampling by the PLANET and GMAN 
collaborations (Albrow et al. 1996, 1997, 1998 and Pratt et al. 1996, Alcock et al. 1997b) 
as well as spectroscopic follow-up of particular events (Lennon et al. 1997). Over 60 events 
are currently alerted per year towards the Galactic Bulge. Since only a handful of these 
are on-going at any given time, monitoring teams can sample events very densely and with 
high photometric accuracy, enabling the detection of many second order effects, including 
-in principle- planetary anomalies. No clear planetary detections have yet been made in 
this way, but preliminary estimates of detection efficiencies show that PLANET, over the 
next two observing seasons, should be sensitive to planetary anomalies caused by Jovian 
planets orbiting a few AU from their parent star (Albrow et al. 1998). Thus, if these kinds 
of planets are common, they should be detected soon. If not, microlensing will be able to 
place interesting upper limits on the frequency of such systems. 

These observational developments have been accompanied by an explosion of 
theoretical work, including further studies of detection probabilities and observing 
strategies incorporating a variety of new effects (Bolatto & Falco 1994, Bennett & Rhie 
1996, Peale 1997, Sackett 1997, DiStefano & Scalzo 1998a,b), demonstration of planetary 
microlensing light curves (Wambsganss 1997), explorations of the degeneracies in the 
fits of planetary events (Gaudi & Gould 1997, Gaudi 1998), and a study of the relation 
between binary and planetary lenses (Dominik 1998). It would thus seem that the 
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theoretical understanding of the detection and characterization of planetary systems using 
microlensing should be well in hand. 

Surprisingly, however, the field still has surprises to offer. Recently, Griest & Safizadeh 
(1998, hereafter GS) came to a rather startling conclusion: for microlensing events with 
minimum impact parameter u m j n ~0.1 (maximum magnification A~10), the detection 
probability is nearly 100% for Jovian mass planets with projected separations lying within 
0.6-1.6 of the Einstein ring of the primary, i.e., the so-called "standard lensing zone." In 
fact, GS found that the detection probability for this subset of events is higher for all 
projected separations, and preferentially so for smaller separations. Since the probability 
that an event will have impact parameter u m i n < 0.1 is ~10%, this means that, for ~10% 
of all events, the existence of a planet in the lensing zone can be detected or ruled out. The 
primary point of this Letter is to stress that the conclusions of GS necessarily imply that, 
regardless of orientation, all Jovian mass planets in the lensing zone dramatically affect 
the central region of the magnification pattern, and thus have a significant probability of 
being detected (or ruled out) in small impact parameter (high magnification) events. 

We present here a preliminary exploration of microlensing by lenses orbited by 
multiple (two) planets. Because our results are intimately tied to those of GS, we refer the 
reader to that paper for a more thorough investigation of detecting single planets with high 
magnification events. Note that we will use high magnification here to mean events for 
which the minimum impact parameter from the primary is u m i n < 0.1. In order to assess 
the frequency with which multiple planets may lie at detectable separations, we calculate 
in §2 the probability of two planets having projected separations in the standard lensing 
zone, and indicate why an even larger zone is more appropriate for high magnification 
events. In §3, we briefly review the formalism needed for calculating the magnification 
patterns created by single, double, and triple lenses, and in §4, we present sample light 
curves. In §5, we discuss possible implications of our results and conclude. 



The "standard lensing zone" is generally defined as the annular region in the source 
plane with 0.6 < r < 1.6Re, where Re is the Einstein ring of the parent star, 



and M is the mass of the primary lens. For the scaling relation on the far right-hand-side, 
we have assumed a source distance Dos = 8kpc and the lens distance Dol = 6kpc. With 
these assumed distances, the lensing zone corresponds to 2.1 — 5.6 AU for a 1.0M Q primary, 



2. 



Lensing Zone 



Frequencies for Multiple Planets 




(2.1) 
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and 1.1 — 3.0 AU for a O.3M primary. As demonstrated by Gould & Loeb (1992), this 
zone roughly corresponds to the range of projected separations b = r / Re for which planets 
will have substantial detection probabilities, averaging over all possible events. However, as 
we will discuss, for the subset of high magnification events, the relevant zone of planetary 
separations is somewhat more extended. 

The standard lensing zone boundaries are defined by the image positions for single 
lens microlensing for a source position equal to Re, the largest source position for which 
current microlensing experiments will alert (magnification A = 1.34). Source positions 
closer to the lens will result in larger magnifications and image positions closer to Re- 
If planetary detection is defined as the source crossing the caustic structure (induced by 
the binary lens) that lies near the planet position, then the planet must be close to these 
image positions and thus within the lensing zone in order to have a measurable effect. 
With such a definition of planetary detection, one would expect multiple planet detection 
to happen only rarely since the source trajectory must traverse both planetary positions, 
both of which must lie within the lensing zone. 

For high magnification events (A > 10), GS have shown that the planets with mass 
ratio g~ 0.001 may be detected with nearly 100% probability even when they lie somewhat 
outside the lensing zone. This is because the planetary anomaly is caused by the source 
approaching or crossing the central caustic (near the primary), not the planetary caustic. 
This in turn implies that the detection probabilities for multiple planets in the lensing zone 
will also be high, providing such a scenario occurs frequently enough. 

Given two planets with true separations (in units of Re) of a\ and a-i, we thus wish 
to calculate the probability that the projected separations b\ and 62 will fall in the lensing 
zone. The relation between the true and projected separations is, 

bk = Ofe[cos 2 4>k + sin 2 fa cos 2 i] 1 ^ 2 , (2.2) 

where i is the orbital inclination, fa is the orbital phase of planet k, and we have assumed 
circular, co-planar orbits. The calculation of the probability involves a three-dimensional 
integral over cosi, fa, and fa, the distributions of which are flat. The result is shown in 
Fig. 1 as a function of 02, for several discrete values of a\ representing known or plausible 
planetary systems with Jovian mass planets. The separations in physical units (AU) scale 
according to Eq. 2.1. 

It is apparent from Fig. 1 that the probability of two Jovian planets falling in the 
lensing zone, regardless of their relative positions, may be quite high. Note, in particular, 
the long tail for higher true separations 02- Furthermore, the conditional probability (lower 
panel of Fig. 1), i.e., the probability that both planets will be in the lensing zone given 
that one of the planets already meets this criterion, is even higher. For high magnification 
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events, this implies that it is highly probable that if deviations from one planet are present, 
deviations from the second planet are present as well. For planets with true separations 
equal to that of Jupiter and Saturn (5.2 and 9.5 AU, respectively), the probability of both 
planets being in the lensing zone is 14% if the planets are orbiting a 1.0M© primary, and 
1% for a 0.3M Q primary. 

Radial velocity techniques have discovered several Jovian-mass planets, many with 
orbital separations substantially smaller than 1 AU (Mayor & Queloz 1995, Butler & 
Marcy 1996), making them difficult to detect via microlensing. Other planets detected 
by radial velocity methods, like the 3Mj mass planet orbiting the GO star 47UMa on a 
circular orbit at 2.1 AU, would be detectable in high magnification events, especially if in 
combination with other planets. As the upper panel of Fig. 1 shows, the planet orbiting 
47UMa would almost never fall in the standard lensing zone of a l.OM primary, but would 
have a probability < 50% of falling within a slightly extended zone defined by 0.5-2.0i?E 
simultaneously with other planets orbiting over a wide range (middle panel of Fig. 1). This 
distinction is important since, as GS have shown, the probability of Jovian-mass detection 
in high magnification events remains high even in this expanded zone. 

The frequency with which multiple planets will reveal themselves in high amplification 
events depends of course on their actual frequency and the distribution of their orbital 
radius and mass. Consider a familiar system of a Jupiter and Saturn orbiting a I.OMq 
primary. Convolving the detection efficiencies of GS as a function of projected separation 
b with the likelihood that Jupiter would have that b simultaneously with Saturn falling 
in the extended 0.5 - 2.QRe lensing zone, we find that ~ 12% of events with u m i n < 0.1 
would reveal the existence of the multiple planets. Since events with u m i n < 0.1 constitute 
~ 10% of all detected events, intense monitoring of 100 events per year could be expected 
to yield ~ 1 multiple-planet lensing event per year, if all Galactic lenses have planetary 
systems like our own Solar System. 

3. Single, Double and Triple Lenses 

In this section, we briefly review and apply the formalism needed for calculating 
the magnification resulting from single, double and triple lens configurations. (For a 
more comprehensive review, see Schneider & Weiss 1986, Paczyhski 1996, and references 
therein.) 

Consider a source with projected position (£,77). Following Witt (1990), we write this 
in complex coordinates as ( = £ + irj. Lensing is the mapping from the source position ( 
to the image positions z = x + iy given by the lens equation, which for N point masses is 
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(Witt 1990): 

N 



k z k -z 



where is the (complex) coordinate of mass k, e& is the fractional mass of lens component 
k, and all distances are in units of the angular Einstein ring, 



9e = 



4GM D LS 
Dql Dqs c 2 



l / 2 R F 

- it, (3 ' 2) 



The magnification Aj of each image j is given by the determinant of the Jacobian of the 
mapping (3.1), evaluated at that image position, 



3 |detJ| 



det^l-||. (3.3) 



In microlensing, the images are unresolved and the total magnification is given by the sum 
of the individual magnifications, A = The set of source positions for which the 

magnification is formally infinite, given by the condition det J = 0, defines a set of closed 
curves called caustics. For the remainder of the discussion, we will label the most massive 
(or only) component of the lens as 1 and define the origin as z\ = 0. 

For the single lens (N = 1) case, the positions and magnifications of the two resulting 
images can be found analytically and their total magnification is 

a u 2 + 2 

u(u z + 4) i / / 

where u = \(\. For u — > 0, Aq — > oo, and the point u = defines the caustic in the single 
lens case. For rectilinear motion, u = \{t — to) 2 /t\ + Umin] 1 ^ 2 , where to is t ne time of 
maximum magnification, u m i n is the minimum impact parameter, and ig; is the time scale 
of the event, 

where v± is the transverse velocity of the lens relative to the observer-source line-of-sight. 
A single lens light curve is then given by F = AqFq, where Fq is the unlensed flux of the 
source, and is a function of four parameters: to, tE, u m i n , and Fq. 

For a double lens (N = 2), Eq. (3.1) is equivalent to a fifth-order complex polynomial 
in z. The solution yields three or five images, with the number of images changing by two 
as the source crosses a caustic. A binary lens generates one, two, or three caustic curves, 
in all cases separate and non-intersecting. The light curve of a binary lens is a function of 
seven parameters: the four parameters describing the single lens case, with the additional 
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parameters b, the separation of the lenses in units of 9e , q, the mass ratio of the system, 
and 9, the angle of the source trajectory with respect to the binary axis. 

Adding a third component to the lens system increases the complexity enormously, in 
particular the caustics can exhibit self- intersection and nesting. Eq. (3.1) is now equivalent 
to a (rather cumbersome) lOth-order polynomial in z. There are thus a maximum of 
ten images, and a minimum of four images, with the number of images changing by a 
multiple of two (Rhie 1997) as the source crosses a caustic. Although in principle the 
image positions can be found by solving the lOth-order equation, the process is slow and 
one cannot easily account for finite source sizes. We therefore adopt here the alternative 
approach of inverse ray-shooting. We first sample the image plane very densely and then 
bin the source positions (Eq. 3.1) in the source plane. The ratio of the resulting density in 
the source plane to the density in the image plane gives the magnification, and repeating 
for all source positions generates a magnification map in the source plane. We then 
convolve this map with different source profiles to produce magnifications appropriate for 
finite sources (Wambsganss 1996). Linear interpolation of the final map yields light curves 
for a particular source trajectory. A triple lens light curve is a function of ten parameters: 
the four single lens parameters, the separations and mass ratios b±, 62 and gi, q 2 , the angle 
of the source trajectory 9, and A9, the angle between the position vectors of the two 
companions. 

For binary and triple systems with small mass ratio(s), most source positions have 
magnifications that are nearly identical to that of a single lens, Aq. It is thus useful to 
define the fractional deviation, 5 = (A — Aq)/Aq, where A is the true (binary or ternary 
lens) magnification. 

4. Illustrating the Effect of Multiple Planets 

An exhaustive study of triple lenses, which would necessitate a exploration of the q±, 
(/2) b±, 62 j and A9 parameter space, is quite beyond the scope of this Letter. However, in 
order to illustrate the effect that a third lens would have on typical light curves we consider 
Jovian planets orbiting stars common in the Galaxy. Fixing b\ = 1.2 and q\ = 0.003, 
corresponding to a Mj planet orbiting a O.3M primary, or a 3Mj planet orbiting a M Q 
primary, we vary the 62 of the second planet with mass ratio q2 = 0.001, corresponding to 
a Saturn-mass planet (for the 0.3Mq primary) or a Jupiter-mass planet (I.OMq primary). 
We concentrate on only those source positions |£ < 0.1 for which the planets have a 
significant cooperative effect. The panels of Fig. 2 show the magnification pattern for 
separations of b 2 = 0.8, 1.0, 1.2, and 1.4 and relative angles A9 = 0,60°, 120°, and 180°. 
For comparison, we also show the magnification pattern when only the planet of mass 
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ratio q± = 0.003 is present. For these maps, we have adopted a uniformly-bright source 
with radius appropriate to a main-sequence star, p = 6*/6e — 0.003(i?*/-R©), where 9* and 
i?* are the angular and physical sizes of the source, and we have assumed Dol = 6kpc, 
D OS = 8kpc, and M = O.3M . 

Note that the case with 62 = 1-2 and AO = is completely degenerate with those 
from a single planet of mass ratio q = qi + q2 = 0.004. While for other configurations the 
magnification patterns with and without the second planet appear dramatically different, it 
should be kept in mind that what one actually measures are light curves, one-dimensional 
cuts through these diagrams. Light curves are shown in Fig. 3, with source radii of 
p = 0.003 and p = 0.01, for the sample source trajectory indicated in Fig. 2. Some 
geometries give rise to light curves that deviate dramatically from the case with only one 
planet, but those with A6 = have shapes that are very similar to single-planet lensing, 
though with larger amplitude and duration. In other words, some geometries with two 
planets (i.e., those with Ad ~ or ~ 180°) will give rise to light curves that are degenerate 
with single planets of larger mass ratios. Furthermore, note that regions of positive and 
negative deviations are more closely spaced when two planets are present. When finite 
source effects are considered, the overall amplitude of the multiple planet anomaly will 
thus be suppressed. Examples can be seen in the 62 = 1-2 and 0.8 and A6 = 180° panels 
of Fig. 3, where for source radius p = 0.01 the amplitude of the anomaly is smaller for the 
double planetary system than the single-planet system, while for p = 0.003 the amplitudes 
are similar. Overall detection probabilities may thus be lower for high magnification events 
when multiple planets and large sources are considered. 

5. Implications and Conclusion 

In this Letter, we have demonstrated that: (1) the probability of two planets having 
projected separations that fall in the "standard lensing zone" (0.6 < b < 1.6) is quite high, 
~ 1 — 15% for planets with true separations corresponding to Jupiter and Saturn orbiting 
stars of typical mass; (2) the influence of multiple planets in and somewhat beyond the 
standard lensing zone can be profound for high magnification events (u m i n < 0.1,) however 
(3) for some geometries, the magnification pattern and resulting light curves from multiple 
planets are qualitatively degenerate with those from single-planet lensing, and (4) for high 
magnification events, finite source effects are likely to suppress more substantially the 
amplitude of multiple planet deviations than single planet deviations. 

Given these results, it would appear that the effects of multiple planets on the 
detection and characterization of planetary systems warrant future study. All previous 
theoretical studies have calculated microlensing planet detection sensitivities either by 
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ignoring multiple planets or by treating each planet independently. For high impact 
parameter events (low magnification), this is probably a fair assumption, but as the 
magnification maps in Fig. 2 illustrate detection probabilities will need to be revised for 
small impact parameters (large magnification). The sense of revision will likely depend on 
finite source effects. It is also likely that for some geometries serious degeneracies exist 
between light curves arising from multiple and single planet high magnification events; 
these degeneracies are above and beyond those present in the single planet case discussed 
by Griest & Safizadeh (1998). This possible degeneracy is especially pertinent in light 
of the fact that the conditional probability of having two planets in the lensing zone is 
substantial. Thus, the interpretation of any given high magnification event may be difficult: 
the degeneracies should be characterized and their severity determined in order to have 
a clear understanding of the kinds of systems whose parameters can be unambiguously 
determined from the deviations. Finally, the calculation of planet detection efficiencies for 
high magnification events should consider multiple planets in order to be able to reliably 
convert the observed frequency of planetary deviations into a true frequency of planetary 
systems. 

We would like to thank the members of the PLANET collaboration, and especially 
Martin Dominik, for comments on an earlier version of this Letter. This work was 
supported by grant AST-95-3061 from the NSF, the Kapteyn Astronomical Institute and 
NWO grant 781-76-018. 
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Fig. 1. — Upper Panel: The probability that two planets with true separations a\ and 02 
(in units of the Einstein ring radius Re of the system) will simultaneously have projected 
separations, b\ and 62, lying in the standard "lensing zone," defined as 0.6 < b < 1.6. 
The probability is shown as function of 02, for a± = 1.5 (solid, appropriate to Jupiter 
orbiting a l.OM primary), a\ = 0.6 (dotted, appropriate to 47UMa) and a\ = 2.7 (dashed, 
appropriate to Jupiter orbiting a 0.3Mq primary). The probability for two planets with 
true separations of Jupiter and Saturn are indicated for a solar mass primary (star) and a 
0.3M© primary (dot). Middle Panel: Same as upper panel, but for the extended "lensing 
zone," 0.5 < b < 2.0. Lower Panel: The conditional probability that both b\ and 62 have 
projected separations in the extended "lensing zone," given that either 61 or 62 satisfies this 
criterion. 
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-0.1 0.1 

Fig. 2. — Contours of constant fractional deviation 5 from the single mass lens magnification, 
as a function of source position (£, rj) in units of the angular Einstein ring, Oe- The 
parameters of planet 1 are held fixed at q\ = 0.003, &i = 1.2, while the projected separation 
62 and the angle between the axes, AO, are varied for a second planet with (72 = 0.001. The 
offset panel is the case when only planet 1 is present. Contours are 5 = ±5% (light lines) 
and ±20% (bold lines). Positive contours are red, negative contours are blue. The caustics 
(5 = 00) are shown in black. A trajectory with minimum impact parameter u m i n = 0.025 
and angle 9 = 260° with respect to axis 1 is shown. 




Fig. 3. — The fractional deviation S from a single mass lens as a function of time for the 
trajectories shown in Fig. 2. The black line is for a source of radius p = 0.003 in units of 
Be', the red line is for p = 0.01. 



